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Abstract 

Dynamics of two-level atoms interacting with their own radiation field in a single-mode 
high-quality resonator is considered. The dynamical system consists of two second-order 
differental equations, one for the atomic SU( 2) dynamical-group parameter and another 
for the field strength. With the help of the maximal Lyapunov exponent for this set we 
investigate numerically transitions from regularity to deterministic quantum chaos in such 
a simple model. Increasing the collective coupling constant b = SnNod 2 /hu we observed 
for initially unexcited atoms usual sharp transition to chaos at b c ^ 1. If we take the 
dimensionless individual Rabi frequency a = 17/2 to as a control parameter, then a sequence 
of order-to-chaos transitions has been observed starting with the critical value a c ~ 0.25 at 
the same initial conditions. 


1 Introduction 

When studying field-matter interactions it is usually of interest to consider the possibility of 
controlling the temporal behavior of the field and/or the atomic subsystems. Say, in resonator 
quantum electrodynamics, it is important to drive the interaction between atoms, moving through 
a cavity, and a quantizied field mode in such a way to attain specified states of the electromagnetic 
field (Fock, coherent, squeezed, and so on) in the cavity and/or desirable states of atoms leaving 
the cavity. It is inexplicitely supposed that we are able, in principle, to attain any desirable state 
(which is accessible, of course, in quantum mechanics) under an appropriate control. 

However, it has been shown in resent years [1] that beyond the rotating- wave approximation 
even the simple model, consisting of N two-level atoms interacting with their own radiation 
field, may demonstrate unpredictable temporal behavior in the sense of deterministic chaos. Our 
previous results [2] have shown that even slight modification in a model, describing this interaction, 
could create dramatic artifacts. The purpose of this work is to treat the routes to deterministic 
chaos in the framework of the dynamical-symmetry approach which has been proved to be useful 
in investigating regular dynamics of a variety of quantum models [3], [4]. 
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2 Dynamical SU( 2) model 

We consider an ensemble of N identical two-level atoms placed in a single- mode high-cjuality 
resonator with the volume V. Each two-level system is described by the SU( 2) Hamiltonian 


H = huRo - hSle(t){R+ + R-), (1) 

in which the operators satisfy the usual commutation relations 

[Ho, H±] = ±H±, [R + , R-] = 2 Ho, (2) 

and u) is the atomic transition frequency that coincides with the resonator frequency. The indi- 
vidual Rabi frequency 0, is given by 

n = dJ T’ (3) 

where d is the dipole moment of the atomic transition. Atoms interact self-consistently by dipole 
interaction with an electric field, whose strength is written in the form 

E(t) = E 0 e(t), (4) 


where E 0 is the constant amplitude and e(t ) the dimensionless variable, 0 < e < 1. 

We treat the field ab initio semiclassically, assuming that it satisfies the usual Maxwell equation 

- + u 2 E - 4 wu> 2 V, (5) 

at 1 

where V = N 0 d < R+ + R- > is the polarisation created by atoms, N 0 = N/V is the density of 
atoms in the resonator. Substituting r = u>t, we can write the eq.5 in the dimensionless form with 
the derivative with respect to r 

£ = -jy < R+ + R- > ■ (6) 


We have introduced following to [1] the dimensionless constant 


6 = 


SnNod 2 
ku ’ 


( 7 ) 


characterizing the energy exchange between the atomic ensemble and the field. 
In addition another dimensionless constant 


a = 


2c 0 


( 8 ) 


will be used to investigate transitions from order to chaos in our model. The expression (8) is 
simply the dimensionless individual Rabi frequency. 

In the dynamical-symmetry approach, each two-level atom is governed by the following single 
equation for the 51/(2) complex- valued group parameter [5] 


9 - 



9 +(2ae) 2 g = 0 , </( 0 ) = 1 , 


9 ( 0 ) = 0 . 


( 9 ) 
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The derivatives in (10) are also defined with respect to r. 

Thus we have two coupled oscilators (9) and (6) describing the self-consistent interaction 
between two-level atoms and a single-mode classical field. Rewriting (9) and (6) in the equivalent 
first-order form, we obtain the following nonlinear dynamical system 


X\ = 

2aey 2 , 

x 2 = 

-2aey u 

2/i = 

— yi + 2 aex 2 

! h = 

y\ — 2aexi, 

£ = 

-V 


V = £ =F £(*13/1 + * 22 / 2 )- 

Signs — and + in the last equation of (10) refer to the initially unexcited and excited atoms, 
respectively. 

The atom-field system (10) obeys two conservation laws 

*l + xl + y 2 i+yl = f, (li) 

± J~i( x 1 + x l - v\ - vl) ~ (c 2 + V 1 ) ± ^e{x x y x + x 2 y 2 ) = const. (12) 

It should be noted that the variables x x = Reg and x-i = Img are not independent [5]. Therefore 
we have three independent real variables, that is the minimum required for chaos [6], 

For two-level atoms the dynamical system (10) is equivalent to the usually adopted Maxwell- 
Bloch equations. Let us introduce the components of the Bloch vector 


u{t) 

— c \ c 2 T c l c 2’ 


v(t) 

= t(cjc 3 -c,c3), 

(13) 

w{t) 

= |c 2 | 2 -| Cl | 2 , 


where c\ and c<i are the probability amplitudes of lower and upper states respectively. 0 
other hand these components can be expressed in terms of the variables x and y as follows 

n the 

u — 

2(xij/i + x 2 y 2 ), 


V — 

2(xij/2 - X 2 y x ), 

(14) 

W = i/1 

2 -|- y 2 2 - X! 2 - X 2 2 . 


Thus we can rewrite (10) in the standard Maxwell-Bloch form 


u 

= -V, 


V 

= u + 4 a£w. 


w 

= —4 a£v. 

(15) 

£ 

= - v , 

V 

= e~h u - 
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3 Numerical results 


Our model possesses two control parameters a and b. We will numerically treat here transitions 
from order to chaos varying one of them in a certain range and keeping another constant. Chaos 
will be diagnosed with the help of the maximal Lyapunov exponent A, which is a quantitative 
characteristic of deterministic chaos describing the mean exponental rate of divergence of two 
initially adjacent trajectories in a phase space [6]. The sign of A gives up a reliable criterion to 
distinguish between regular and chaotic dynamics of a system in question. When it is neglibly 
small the motion is said to be regular. If A becomes positive for a certain range of values of a 
control parameter a system is chaotic for this range. Chaos may also be confirmed by continuous 
power spectra. 



b 

Fig. 1 The maximal Lyapunov exponent as a function 
of the control parameter b for initially unexcited 
atoms, a=0.25 ■ 1CT® . 
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Fig.2 The maximal Lyapunov exponent as a function 
of the control parameter b for initially excited 
atoms, a=0.25 ■ 10* 6 . 




Fig.3 The maximal Lyapunov exponent as a function 
of the control parameter a for initially unexcited 
atoms, b=1 . 


Fig.4 The maximal Lyapunov exponent as a function 
of the control parameter a for initially excited 
atoms, b=1. 
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By varying the collective coupling parameter b we, in fact, change the density oi atoms A^o in 
a cavity. Numerical integration shows that the maximal Lyapunov exponent A becomes positive 
when b exceeds a critical value b c . Its magnitude depends essentially on initial conditions. It is 
seen from Fig.l that b c ~ 1 for initially unexcited atoms. For initially excited atoms (Fig.2) the 
maximal Lyapunov exponent becomes positive for much smaller critical value of b. 

We have observed a quite different transition to chaos with chaotic regimes alternating among 
regular regimes when varying the individual dimensionless Rabi frequency a and fixing the pa- 
rameter b. Fig. 3 and Fig. 4 demonstrate such a behavior for initially unexcited and excited atoms, 
respectively. 


4 Outlooks 

We have demonstrated two possible routes to chaos in the interaction of two-level atoms with their 
own radiation field. From a more general points of view, we have observed numerically order-to- 
chaos transitions in the system of two coupled nonlinear oscillators (6) and (9). At last, from an 
abstract point of view, we have treated such transitions in a system consisting of the "driven” 
su( 2) group treated as a nonlinear dynamical system. Thus, the results, obtained in this work, 
are applicable in a more general context. They may be applied with slight modifications to any 
driven physical (chemical, biological, ecological, etc) system with the underlying SU(2) dynamical 
symmetry. 
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